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Abstract. We present a supersymmetric analysis for the oi-dimensional Schrodinger equation with the 
generalized isotonic nonlinear-oscillator potential V(r) = B 2 /r 2 + ui 2 r 2 + 2g(r 2 — a 2 )/(r 2 + a 2 ) 2 , B > 0. We 
show that the eigenequation for this potential is exactly solvable provided g = 2 and (loo 2 ) 2 = B 2 + (I + (d — 
2)/2) 2 . Under these conditions, we obtain explicit formulae for all the energies and normalized bound-state 
wavcfunctions. 
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1. Introduction 

Recently, Carinena et al. pQ studied Schrodinger's equation with a quantum nonlinear oscillator potential of the 
form 



d 2 , 2a; 2 - 1 
1 x 2 1 " 



ip n (x) = 2E n ip n (x) (4) 



dx 2 ■ ~ ■ - (2aa + 1)aJ ^(*)=^(*)- (1) 

Their interest in this study came from the fact that Eq. ([]]), is exactly solvable, in a sense that the eigenenergies 
and eigenfunctions can be obtained explicitly; they were able to show that [T] 

(2) 

E n = -3 + 2n, n = 0,3,4,5,..., 
where the polynomial factors P n (x) are related to the Hermite polynomials by 

f 1 ifn = 

Pn(x) = I (3) 

[ H n (x) +AnH n - 2 {x) +4n(n- 3)H n -i(x) if n = 3, 4, 5, . . . 

Soon afterwards, Fellows and Smith 2 considered another interesting case and, in particular, they showed that 
for certain values of the parameters w and g, the potential in the Schrodinger equation 

d 2 2 *^ ^ 

dx 2 (x 2 + a 2 ) 2 

is indeed a supersymmetric partner of the harmonic oscillator potential. By means of the supersymmetric 
approach, Fellows and Smith [5] were able to construct an infinite set of exactly soluble partner potentials, 
along with their eigenfunctions and eigenvalues. Very recently, Sesma [3], using a Mobius transformation, was 
able to transform Eq.(j4]) into a confluent Heun equation and thereby obtain a simple and efficient algorithm to 
solve the Schrodinger equation Q numerically, no matter what values are chosen for the parameters w and g. 
Furthermore, using suitable mass distributions, the position-dependent effective mass Schrodinger equation has 
been solved for a new nonlinear oscillator [3]. We note that the term 2(x 2 — a 2 )/(x 2 + a 2 ) 2 in Eq.Q can be 
written as the sum of two centripetal barriers in the complex plane 

2{x 2 -a 2 ^ 1 i 



(x 2 + a 2 ) 2 (x + ia) 2 (x — ia) 21 

that is to say, a rational potential with two imaginary poles symmetrically placed with respect to the origin [1] . 

The purpose of the present work is to use the supersymmetric approach to analyze the exact analytic solutions 
for the more general potential class We consider the d— dimensional Schrodinger eigenvalue problem 

HV(r) = £*(r), H = -A + V, (5) 

where A is the d— dimensional Laplacian operator, d > 2, and V(r) is the central potential given by 

V{r) = ^+u 2 r 2 +2g f- a l ) B 2 >0. (6) 

We show, through the factorization method, that the eigenequation for this potential, Eq. ([S]) , is indeed exactly 
solvable provided g = 2 and (wo 2 ) 2 = B 2 + (£ + (d — 2)/2) 2 . The article is organized as follows. In section 

2, we consider the d-dimensional Schrodinger equation with the nonlinear oscillator potential (j6]). In section 

3, we briefly review the factorization method and the basic fomulae from supersymmetric theory that we need 
for our investigation 5,6 . In section 4, we show that the isotonic nonlinear-oscillator potential p][2| and 
the Gol'dman and Krivchenkov potential are isospectral supersymmetric partner potentials. In section 6, the 
explicit construction of all the exact solutions of Schrodinger's equation ([5]) with V(r) given by © is presented. 
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2. Schrodinger's equation in d dimensions 

By considering the action of the Laplacian operator on a wavefunction of the form ty(r) = u(r)Yi(9i,02, ■ • ■ , Od-i) 
with a spherically symmetric factor u(r) and a generalized spherical harmonic factor Yi, we obtain the radial 
Schrodinger equation in d dimensions for a spherically symmetric potential V(r) as 015]: 

d 2 u d-1 du l{l + d-2) r 
dr 2 r dr r 2 
where, if r € 3? d , then r = ||r||, and 

u{r) e L 2 ([0,oo), r d_1 dr). 
The first-order derivative term can be removed by using the new radial function ip(r) given by 



-it + V(r)u = Eu, 



(7) 



Thus, we find 



where 



^(r) =r-a-«(r), d > 2, ^(0) = 0, ?^(r) 6 L 2 ([0,oo),dr). 

d 2 V> 



dr 2 



+ [/(r)V' = £ty 



U(r) = 7(r) + 



(21 + d- l)(2l+rf-3) 
4^ ■ 



We now recall the definition of V(r) in and combine the terms in 1/r 2 to obtain, 



U(r) 



k(k + 1) 



+ uj 2 r 2 + 2g 



(r 2 - a 2 ) 



where k is defined so that 
1 



B 2 



(r 2 + a 2 ) 2 ' 
(21 + d - l)(2i + d-3) =k(k + l), 



(8) 
(9) 
(10) 



that is to say, 



B 2 



d-2 



(11) 



We shall find it convenient to label energy eigenstates ip n k = ipw/d, where n = 0, 1, 2, ... is the number of radial 
nodes, and 7^ = k + h. 



3. The factorization method 

In this section we give a brief review of some concepts of supersymmetric quantum mechanics (SUSY QM) 
that we shall need in the following sections, namely the factorization method, and supersymmetric partner 
potentials. We start with Schrodinger's time-independent equation for a one-dimensional radial problem (in the 
units h — 2m = 1) 



Hip(r) = 



tjj(r) = Eip(r) 



(12) 



where the potential U(r) is real and possibly singular at r = 0. The wave function ip( r ) must be square-integrable 
and normalized in a sense that J \tp(r)\ dr = 1. The main idea of the factorization method, as introduced in this 
context by Schrodinger and Dirac, is to write the second order differential operator H in (fT2j) as the product of 
two first order differential operators and A, such that H = A^A and 



A f = -4- +W(r), A = ^-+W(r) 
dr dr 



(13) 
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where the function W(r) is a real function of r and is known as the superpotential of the problem. The operators 
A and are Hermitian conjugates of each other: (A^)* 1 = A. Let us define 

Hi = -4U = ~ + EMr), #2 = AA* = + W (14) 

In general Hi = H and H 2 are two different Hamiltonians. They are known as partner Hamiltonians in SUSY 
QM and are given explicitly by: 



and 



IT, . JtA = (-| + WD) (I + H.- W ) . -£ + W\r) - *™ - -£ + t«r) (15) 

* - = (s + B ' M ) (-s + H,( ->) ~£ + + -£ + <16) 



where the potentials C/i and E/2 are known as supersymmetric partner potentials defined by: 



t/i(r) = IU» and f/ 2 (r) = W 2 {r) + — (17) 

dr dr 

Suppose that ipn\r) and ipn\r) are the eigenfunctions of Hi and _ff 2 , respectively. In the unbroken SUSY 
case, the ground state is not degenerate with a vanishing energy E = and it is usually expressed in terms of 
the superpotential W(r) as: 

A 4 1] (r) = (J: + W ( r )) 4 1} (r) = => ^ (r) = C cxp | W(r)dr^j , (18) 

where C is the normalization constant. The key result is the iso-spectrality between the two Hamiltonians H\ 
and H 2 for all but the ground state (n = 0), which can be shown as follows. Since the energy eigenvalues of 
Hi and H 2 are positive semi-definite E^ 2 > 0, we have for H 2 i\ffi — AA^ipn^ = En^ipn^ and by multiplying 
through by A^ we see that 

HiiAHP) = AU(AW*) = Ah 2 1><? = E^(AH^). (19) 

Thus A^tpn^ is an eigenstate of Hi with same energy eigenvalue En \ and there must be ipn\r) such that 

^(r) = c n A^i 2 \r) (20) 

where c„ are constants for n — 1,2,... Similarly, for HiipiP = A^AipiP = En^ipn^ and, multiplying through 
by A, we find 

H 2 (A^P) = AAHAtPP) = AHrfW = E£\A^). (21) 
Thus AipiP is an eigenstate of H 2 with the same eigenvalue En \ and there must be ipn\r) such that 

^\r) = c' n A^\r). (22) 
Furthermore, since 

1 = (V^WIV^W) - \c' n \ 2 {A^\r)\A^\r)) = \c' n \ 2 {ip^ (r) \A f Aip^P (r)) = ^I^MIV^M), 



we have 



^)(r) = (^-^(r) = (Elpyi/ 2 (± + TU(r)) ^«(r). (23) 
Similarly, we find 

^)(r) = (£( 2 ))- 1 /Mt^ 2 )(r) = (Ei?y 1/2 (-^ + W(r)) V^M- (24) 

Thus, if we knew the eigenvalues and eigenfunctions of either of the two partner potentials, we could immediately 
derive the spectrum of the other. However, the above relations only give the relationship between the eigenvalues 
and eigenfunctions of the two partner Hamiltonians, but do not allow us to determine their spectra. In the 
next, we are guided by the idea of finding pairs of (essentially) isospectral Hamiltonians, one of which has a 
known soluble Hamiltonian. 
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4. Supersymmetric partner of the Isotonic potential 

A key step in this work was the idea to consider a candidate for a superpotential W(r) by means of the following 
Ansatz: 

W(r) = -+u/ r+-^— I , (25) 
r r A + or 

where k',ui' and s are real parameters to be determine shortly, and a is a fixed potential parameter. The 
potential U\, then reads 

= ^ + + + 8 (-2 a V + 8 + 2^ + l )ra - MW-2^ + l) + + ^ _ 

r Z _|_ fl 2j2 

(26) 

and the potential t/ 2 is given by 

. . , 22 fc'(jfe'-l) s(s-l + 2/c'-2aV)r 2 +a 2 s(l + 2fc' - 2oV) ,„ 

= -'V + + J (r2 j +fl2)2 " + (2* + 2fc' + l)u/ 

(27) 

If we now compare the potential U\ (r) with 

k(k + l) 22 „ (r 2 -a 2 ) 

r 2 ^2 _|_ fl 2j2 

we have for k' = k, J = u, s(-2a 2 uj' + s + 2k + 1) = 2g, and s(2u/a 2 - 2k' + 1) = 2g, that 
s = 4(uja 2 - k), 

(29) 

g = 2(wa 2 - k)(2wa 2 + 1 - 2k). 

With these values of s and g we may reduce £/ 2 by the assumption that s(s — 1 + 2k — 2a 2 uj) = and 
a 2 s(l + 2fc-2a 2 w) = to 

V 2 (r) = cu 2 r 2 + fc(fc ^ X) + (8(wa 2 - k) + 2k + l)w. (30) 
These assumptions are valid under the following conditions: 

2 1 
toa — k + -, g = 2 and s = 2. 



W = ^+-V + 2,AL^ (28) 



In summary, we have two partner potentials 



. i fc(fc + l) 22 4(r 2 -a 2 ) . . . 

W = + wV + ( ,2 +fl2)2 + ^( 2fc + 3 ) ( 31 ) 

V 2 (r) = ^^+u, 2 r 2 +u,(2fc + 5). (32) 



provided the parameters w and a satisfy 



,2_, , 1 ,/o2 , ( t , (d-2) 



2 



ca 2 = fc + - = W J B 2 +(£+^ r ^) . (33) 
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5. Gol'dman and Krivchenkov potential in d-dimensional 

Schrodinger's equation with the Gol'dman and Krivchenkov potential [9] in d-dimensions is given by 



A(A + 1) 



dr 2 



+ Pr 2 + ^ ) ipn-yAr) = En^ipn-yAr), A = i + -(d-3), d>2, (34) 



and has exact eigenvalues given by [5] 



E, 



2^(2n + ld ), »,Z = 0,1,2,..., 



(35) 



where 



ld = l + \ja + {k+-) 2 . 



Meanwhile, the exact eigenfunctions are given by [5] 



'2/3W2( 7d ) n 



n\T{ ld ) 

Here, we use the Pochhammer symbol ("f) n , where 



r -w-i/2 e -iv7»r ^(-n^V^r 2 ). 



(36) 



(7)n = 7(7+l)(7 + 2)...(7 + n-l) 



r( 7 + n) 

r(7) • 



and i-Fi is the confluent hypergeometric function denned by ( [10] . chapter 7) 

n 

1 F 1 (-n;b;x) = ■ 



(6)* fc! ' 



(37) 



Using these exact solutions for the Gol'dman and Krivchenkov potential, we can show that the Schrodinger 
equation 

d 2 + k(k-l) 
dr 2 r 2 



+ w 2 r 2 + w(2k + 5) ip n k(r) = e nk ip nk (r) 



has exact solutions given by 

e nfe = 2w(2n + 2k + 3), 

and 



, / n , 1W /2w fc + 1 /2(fc + l/2)„ , 12 ^ ^ 1 2N 

^ fc(r) = ( - 1} V n!r (fc + l/2) r e 2 ^(--^+2-- 2 )> 



(38) 
(39) 

(40) 



where 



k + - = 7d = wa 2 = 



B 2 + \l 



d-2 



(41) 



It is clear from these relations that states having the same value for the combination 2£ + d are degenerate. 
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6. Exact solutions for isotonic nonlinear-oscillator potentials 

The supersymmetric approach, briefly discussed in section 3, along with exact solution of the Gol'dman and 
Krivchenkov potential allow us to obtain the exact solutions of Schrodinger's equation 
d 2 k(k + 1) Mr 2 — a 2 ) \ 

j-2 + -^T A +^ 2 + ( , 2+Q2) 2 +^( 2fc + 3 )j Mr) = e nk Mr), (42) 
namely, 

e nk = 2w(2n + 2fc + 3), (43) 
and the corresponding wavefunctions are given by 

Mr) = C' n (-± + t +ujr + _^L\ r k e-i" r2 LtHur 2 ) (44) 



dr r r 2 + a 2 



where 



C = 2^ + i/ 2 (fc + l/2)„ 

1 ; Y 2w(2n + 2fc + 3) n!r(fc + l/2)' { ' 

Here we have used the well-known relation between the confluent hypergeometric function (j37]) and Laguerre 
polynomials ( [10], page 203) 

n 1 

1 F 1 (-n;a+l;z) = — L«{z). 

(a + 1)„ 

A straightforward computation, aided by the differential identity of Laguerre polynomials ( [TO], page 203, 
formula (11)), namely 

±I%{z) = -Kt\{z), (46) 

yields 

(~<l r k+l -ur 2 /2 f 

^(r) = ^Vf^ ( 2k + 2n + 3 ) L "r V W) - 2 (n + V 2 ) + 2uj a 2 L k ^l 2 (ur 2 ) 

r z + a z 

(47) 

where we have used the identity ( [10], page 203, formula (8)) 

L%(z) = L%_ 1 (z) + L%- 1 (z). (48) 
Since wa 2 = k + 1/2, we may write Eq . (j47|) as 

(2k + 2n + ^Ll- 1 ' 2 ^ 2 ) - 2(n + I) 2 {cur 2 ) + {2k + \)L k n +1 / 2 {ur 2 ) 

(49) 



, , , C' n r k+1 e-^ 2 l 2 

<Pnk{r) = 5- o 



This may be compared with the results of Fellow and Smith [5] by use of the well-known relations between 
Hermite and Laguerre polynomials ( |10) . page 216, problem (1)) 

(— 1 \ n 

L " 1/2( " r2) = H2n[ ^ T) (50) 

and 

L)' 2 {Lor 2 )= \ / H 2n+1 (V^r). (51) 

For k — 0, w = 1 using the identity ( [TOJ, page 188, formulas (4) and (6)) 

2zH k (z) = H k+l {z) + 2kH k ^(z). (52) 
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we obtain 



4>nk {r) 



{-\TC n e~r 2 n 



4(2n + 3)J?2n+i(r) + 8n(2n + 3)F 2 „-i(r) + # 2 «+3(r) , 



(53) 



2 2 "+ 2 n! (r 2 + a 2 ) 



where 



C = (-1) 



I) 




(54) 



(2n + 3) n! tt' 



These general results agree with the results discussed by Carihena et al pQ, and Fellows et al. [2] for the 
comparable odd solutions n = 3, 5, 7, ... in Eq.Q. 

7. Conclusion 

In this paper, we have studied a family of generalized quantum nonlinear oscillators in d-dimensions. These 
isotonic oscillator problems have potentials V(r) of the form 



We show that if g — 2 and (wo 2 ) 2 = B 2 + (£+(d— 2)/2) 2 , then this potential can be regarded as a supersymmetric 
partner of the Gol'dman and Krivchenkov potential, for which exact solutions can be constructed. Thus, under 
these conditions, we are able to solve the eigenproblem Htjj = (—A + V)ip = Eijj exactly in d dimensions, and 
provide formulae for all the discrete eigenvalues and corresponding normalized wave functions. We have shown 
that our solutions to the general problem agree with the results reported earlier in Refs. [1] and [2] for the 
comparable odd states n = 3, 5, 7, . . . in d = 1 dimension with B = 0. 
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V(r) = —+ u 2 r 2 + 2g 



r 



(r 2 - a 2 ) 

( r 2 +a 2)2' 



(55) 



